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1.  Introduction 


Enhancenieut  of  the  maneuverability  and  controllability  of  current  and  future  fighter  aircraft  re«|nirefi 
the  ability  to  more  accurately  predict  the  aerodynamics  of  these  vehicles  during  rapid  maneuvers  and  at  tlu' 
extremes  of  their  flight  envelopes.  This  has  led  to  a  particular  interest  in  high  angle  of  attack  and  unsteady 
aerodynamics.  One  of  the  dominant  features  of  the  flow  field  around  a  fighter  aircraft  at  angle  of  attack 
is  the  vortex  flow  above  the  delta  wing.  A  better  understanding  of  this  flow  field  will  lead  to  an  improved 
understanding  of  the  aerodynamics  of  a  fighter  aircraft.  The  vortex  flow  over  a  delta  wing  also  provides  a 
good  model  for  studying  complex,  three-dimensional  vortical  flows.  This  work  is  part  of  a  continuing  efl'urt 
to  improve  predictive  capabilities  for  both  high  angle  of  attack  and  unsteady  aerodynamics. 

The  present  work  has  developed  out  of  a  series  of  experimental  and  computational  invest  igations  of 
the  following  background.  Gad-el-Hak  and  Blackwelder  [1]  observed  for  water  tunnel  tests  of  two  della 
wings  that  the  vortex  sheet  which  emanates  from  the  leading  edge  of  the  delta  wing  rolls  up  into  discrete 
vortices  that  undergo  a  pairing  process.  The  instability  appears  as  alternating  dark  and  light  regions  in 
the  dye  sheet  near  the  leading  edge  of  the  delta  wing  and  extending  along  it,  figure  1.  This  instability  is 
said  to  be  similar  to  the  instability  and  pairing  process  described  by  Brown  and  Roshko  [2]  and  Wiiunit 
and  Browand  [3]  for  two-dimensional  shear  layers.  Both  are  related  to  the  classical  Kelvin-Helmholtr.  type 
instability  of  two-dimensional  shear  layers. 

Payne  d  al  [4]  also  observed  an  instability  in  the  shear  layer  forming  the  primary  vortex  usii'.g  smoke 
flow  visualization  technieptes.  They  again  relate  the  growth  of  these  secondary  structures  to  the  cla.s.sii  al 
Kelvin-  Helmholtz  type  of  instability..  The  v.bserved  vortical  structures  appeared  to  be  static,  howev«*r,  and 
did  not  rotate  with  the  vortex.  They  relate  these  observations  [5]  to  those  of  Gad-el-Hak  and  Blackwelder. 

Lowson  [6]  ijas  carried  out  low  speed  wind  tunnel  tests  to  investigate  the  vortex  sheet  structure  formed 
about  delta  wings.  In  this  work  he  concludes  that  not  one,  but  two  types  of  instability  exist  in  the  shi*ar 
layer  leaving  the  leading  edge.  The  first,  shown  in  figure  2  is  the  quasi  two-dimensional  type  of  instability 
observed  by  Gad-el-Hak.  The  second  is  a  locally  streamwise  instability  of  the  vortex  sheet  that  gives  rise 
to  steady  structures  within  the  vortex  core.  Lowson  suggests  that  this  is  the  instability  obs<?rved  by  Payne 
(I  nl. 

All  tiu'se  experimental  results  have  been  obtained  for  scaled  models  at  relatively  low  Reynolds  number. 
Recent  photographic  evidence  shows  that  similar  types  of  secondary  vortical  structures  are  observed  for 
the  LEX  vortex  flow  on  an  F-18  in  flight  (7j.  This  strongly  indicates  that  the  vortex  structure  found 
experimentally  at  lower  Reynolds  number  occurs  in  fliglit  conditions  at  full  .scale  Reynolds  number. 


1 


A  .similar  tjpe  of  shear-layer  instability  was  observeJ  by  the  present  authors  [8]  in  calculations  for  a  7(5® 
sweep  delta  wing  at  20-5®  angle  of  attack  and  a  Reynolds  number  of  900,000.  An  unsteady  behavior  of  the 
shear  layer  which  emanates  from  the  leading  edge  of  the  delta  wing,  related  to  the  instability  .!escril>ed  by 
(Jad-el  Hak  and  Blackwelder.  was  found  in  these  calculations.  Uncertainties  in  the  effect  of  grid  i-esoliitioii. 
temporal  accuracy,  and  Reynolds  number  were  noted. 

'I'he  current  work  addresses  some  of  the  uncertainties  noted  for  the  previous  work.  Further  coiiipiitatioiis 
are  performed  at  a  lower  Reynolds  number,  Re  =  .50,000.  and  for  enhanced  griti  resolution.  An  improved 
di'scription  of  the  unsteady  process  observed  will  be  given.  The  lower  values  of  Reynolds  number  will  also 
provide  the  ability  to  more  closely  relate  the  minierical  results  to  the  existing  experimental  ob-serval ious 
just  desciibed. 
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2.  Governing  Equations 


The  governing  equations  for  the  present  problem  are  the  unsteady,  three-dimensional ,  full  Navior-Stok«*» 
equivlions  written  in  strong  conservation  form  [9]  using  a  general  coordinate  transformation 
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The  system  of  equations  is  closed  using  the  perfect  gas  law.  Sutherland’s  formula  for  viscosity.  aii<l  i.i.- 
a'sumplion  of  a  constant  Praiidtl  number.  Pr  =  0.72.  Flow  quantities  have  been  nun<limensionalized  oy 
their  respective  free  stream  values  e.\cej>l  for  pressure,  which  is  nondimensionalized  by  twice  (he  freest  ream 
dynamic  prt.>s.sure,  and  speed  of  sound  which  is  nondimensionalized  by  (he  freestream  velocity.  All  leiigllis 
have  been  normalized  by  the  root  chord  length  of  (he  delta  wing. 
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3.  Numerical  Procedure 


TliP  governing  equations  are  solved  numerically  using  the  implicit  approximately  factored  algorilliiji 
of  Beam-Warming  (10).  The  equations  are  differenced  using  Euler  implicit  time-differencing  and  secoinl- 
order  accurate  central  differences  for  all  spatial  derivatives.  A  blend  of  second  and  fourth  order  nonliiu  ar 
dimipation,  as  suggested  by  Jameson  [11],  is  used  to  stabilize  the  central  difference  scht'me.  The  current 
work,  in  which  subsonic  flow.s  are  investigated,  requires  only  fourth-order  dissipation. 

A  fully  vectorized,  time-accurate  solver  has  been  developed  to  implement  the  aforementioned  .sclienu*  (12). 
The  computational  requirements  of  the  scheme  are  approximately  38  word.s/grid  point  and  4.4  x  10“” 
CPU  .sec/grid  point/iteration  on  a  Cray  2.  This  code  has  been  validated  for  a  number  of  simple  laminar 
flows.  Solutions  for  the  case  of  a  supersonic  delta  wing  have  been  compared  with  both  experimental  resiilt.s 
and  other  numerical  results  [13].  Good  agreement  for  both  comparisons  has  been  obtained.  The  code  has 
also  been  shown  to  capture  successfully  transition  from  Couette  flow  to  Taylor- vortex  flow  due  to  centrifu¬ 
gal  instability  [12].  Ibisteady  results  have  been  obtained  for  a  pitching  slender  body  of  revolution  [1 1]. 
Filially,  the  code  has  been  used  to  calculate  both  steady  and  unsteady  horseshoe  vortex  flows  occurring 
at  a  cyliiuler/flat  plate  juncture  [15].  These  results  have  provided  a  good  validation  base  for  the  code.  It 
is  of  particular  importance  that  the  latter  results  demonstrate  the  ability  to  correctly  simulate  naturally 
ocruring  instabilities  in  the  flow  Held  and  the  transition  of  the  flow  from  steady  to  unsteady  behavior. 

Mac('-orinack[16]  suggests  that  a  desirable  form  of  a  numerical  method  is; 

{numeriesjAq  =  {physics}  (3.1 ) 

This  form  allows  modifications  to  the  left-hand  side  of  equation  (3.1)  for  numerical  reasons  (e.g..  efficiency) 
without  affecting  the  accurate  representation  of  the  governing  physical  equations  contained  on  the  riglit.- 
hand  side.  The  current  numerical  procedure  has  this  form  for  steady  flows.  Implementing  a  Newlon- 
like  subiteration  procedure  recovers  this  form  for  unsteady  flows.  This  subiteration  procedure  has  Ijm-h 
successfully  used  by  several  authors  [17.18,19].  The  Newton  subiteration  procedure  can  provide  distinct 
improvements  to  the  current  numerical  scheme.  Subiteration  improves  the  accuracy  of  the  scheme  via 
reduction  or  elimination  of  both  factorization  and  linearization  errors.  Relaxation  of  stability  limits  in  the 
three-factored  scheme  due  to  the  reduction  of  the  factorization  error  improves  the  efficiency  of  the  sclienii . 
Finally,  one  may  obtain  higher-order  temporal  accuracy  in  a  straightforward  fashion. 

'I’he  subiteration  procedure  is  implemented  in  the  currently  existing  code  with  a  few  simple  modinca- 
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tions.  The  standard  Beam- Wanning  algorithm  for  equation  (2.1)  is  written  as  follows; 


=  -AtR’'  (3.2) 

with: 

R"  =  6(F  +  6^G  +  6(H~  +  KG,  +  6^ H,] 

Aq"  =  -  ?'• 

f.(  =  (•{(•■)  + <?i() 

ill  =  K(F  +  Q2,i) 
i<:  =  ^(;(i*  +  Q3<) 

where  .-i.fl.  and  ('  are  the  inviscid  flux  Jarobians  and  (ji.Qs.  and  Qa  are  the  appropriate  visroiis  flux 
Jarobians.  For  the  .subiteration  procedure  one  term  is  added  to  the  right  hand  side  of  equation  (3.2)  nnd 
the  definition  of  Aq"  is  modified,  yielding: 

[j^ + Aa^]r*^[j^ + (-‘S) 

tvhere: 


q  =  q/J 

Aq^  = 

Atf  =z  subiteration  time  step 

Atfhs  =  physical  time  step. 

Here  p  is  the  subiteratioii  count.  For  p  =  0,  qf  =  q'*  and  the  scheme  reverts  to  the  noniterative  form  of  the 
algorithm,  equation  (3.2).  As  p  —  oc,  q^  —  and  the  .solution  at  the  new  time  level  is  obtained. 

CUirreiitly,  the  subiteration  procedure  is  implemented  using  the  same  factorization  .scheme  as  the  non¬ 
iterative  algorithm.  Since  the  physics  of  the  problem  now  is  contained  in  the  right-hand  side  of  i*qua- 
tion  (3.3).  the  left-hand  side  may  be  modified  to  improve  the  efficiency  of  the  numerical  scheme  without 
adversely  affecting  the  time  accuracy  of  the  problem.  This  area  remains  for  future  investigation. 


4.  Grid  Structure  and  Boundary  Conditions 


The  coinputalioiial  grid  for  the  delta  wing  was  obtained  using  simple  algebraic  grid  gener.-?tton  U.*ch- 
iiiqiies.  The  grid  structure  is  of  an  H-H  type  with  planar  grids  being  stacked  in  Mie  axial  direct  ion.  The  H-H 
structure  provides  for  good  resolution  of  the  sharp  leading  edge  with  no  rounding.  All  farfield  boundaries 
are  located  1.5  to  2  chord  lengtlis  away  from  the  delta  wing.  The  grid  configuration  is  .shown  in  figure  -1. 

The  boundary  conditions  are  implemented  in  the  following  manner.  On  the  wing  surface,  the  i.o- 
slip  condition  for  the  velocities,  adiabatic  wall  temperature  condition  and  zero  normal  pressure  gradient 
condition  are  used.  At  the  dowastream  boundary,  flow  variables  are  extrapolated  from  the  interior.  1‘his 
allows  for  convection  of  the  vortices  out  of  the  computational  domain.  .Symmetry  conditions  are  imposed 
along  the  mid-plane  of  the  wing. 

Characteristic  boundary  conditions  [20]  are  used  at  the  upper,  lower,  side,  and  upstream  boundaries. 
Locally  one-dimensional  flow  is  assumed  from  which  the  corresponding  characteristic  variables  are  obtained. 
The  iKiuiidary  conditions  are  then  applied  to  the  characteristic  variables  with  their  values  l)eing  specified 
from  either  the  freestream  or  interior  v-alues  depending  on  the  sign  of  the  corresponding  eigenvalue.  For 
8ul)sonic  inflow  four  characteristic  variables  are  specified  from  the  freestream  conditions  with  the  fifth 
extra|>olateti  from  the  interior.  For  subsonic  outflow  four  characteristic  variables  are  specified  from  the 
inb'rior  svith  one  variable  specified  from  the  freestream  values.  In  the  present  calculations,  the  n|)slream 
and  lower  boundaries  are  inflow  boundaries,  and  the  upper  and  side  boundaries  are  outflow  boundaries. 
The  specification  of  characteristic  conditions  reduce  undesirable  reflections  from  the  boundaries  of  the 
roniputational  domain. 


5.  Results 


5.1  High  Reynolds  Number  Case 

All  results  presented  are  calculated  for  a  76**  sweep  delta  wing  of  aspect  ratio  1  (lig.  3).  This  cunligiira- 
tion  has  been  investigated  experimentally  by  Hummel  [21,22].  The  first  calculations  for  this  configurat  ion  [i<] 
are  made  at  Reynolds  number  900,000,  angle  of  attack  of  20.5‘,  and  a  freestream  Mach  number  0.2.  Th(‘se 
conditions  correspond  to  the  experiment  of  Hummel,  except  for  the  Mach  number  which  is  high(‘r  than 
the  experimental  Mach  number,  0.05.  A  node  system  consisting  of  60  x  151  x  100  grid  points  is  used.  37 
points  in  the  axial  direction  and  80  points  in  the  spanwise  d'l  ‘ction  are  located  on  the  wing.  The  minimiiin 
spacing  at  the  surface  is  A:  =  0.00001  and  a  uniform  spacing  along  the  leading  edge,  Ay  =  O.OOOi  I.  is 
specified. 

The  numerical  calculations  fur  this  case  were  found  to  be  unsteady  unless  a  coarse  grid  is  employed.  'I'lu' 
unsteadiness  in  the  calculations  is  attributed  to  a  Kelvin- Helmholtz  type  instability  of  the  shear  layer  which 
emanates  from  the  leading  edge  of  the  delta  wing  and  forming  the  primary  vortex.  This  unsteady  ficv 
field  is  characterized  by  the  formation  and  shedding  of  small-scale  shear-layer  vortical  structures  from  the 
leading  edge  of  the  delta  wing.  It  is  suggested  that  this  behavior  is  similar  to  that  observed  by  Chul-el-lluk 
and  Blackwelder  in  their  water  tunnel  experiments  [1]. 

A  comparison  of  the  average  surface  pressure  at  several  axial  locations  along  the  wing  with  the  e.\pi'ri- 
mental  data  of  Hummel  [22]  is  given  in  figure  5.  The  calculated  results  compare  well  with  the  experimental 
data  showing  only  a  slightly  greater  expansion  in  the  region  below  the  core  of  the  primary  vortex.  The 
overprediction  of  the  pressure  at  *  =  0.9  is  consistent  with  the  calculations  of  Thomas  [23]  and  is  due  in  part, 
to  the  onset  of  transition  near  the  trailing  edge  in  the  experiment.  Further  details  of  these  computations 
may  be  found  in  reference  [8]. 

An  investigation  of  the  effect  of  time  resolution  on  the  unsteady  behavior  of  the  high  Reynolds  number 
case  is  carried  out.  Two  values  of  At  are  considered,  At  =  0.001  and  At  =  0.00025.  Three  subiteratious 
are  used  for  each  step  at  At  =  0.001.  The  pressure  is  monitored  at  points  traversed  by  the  small  scale 
vortices  shed  from  the  leading  edge.  Results  for  both  time  steps  at  an  axial  location  r  =  0.3  are  given  in 
figure  6.  This  location  provided  the  highest  computed  frequency  and  should  therefore  be  the  most  sensitive 
to  time  step  resolution.  The  frequency  for  the  smaller  time  step,  St  ss  24,  is  higher  than  the  fre(|ueiuy 
for  the  larger  time  step,  St  »  20.  The  basic  nature  of  the  flow  remains  the  same,  however.  A  time  step 
of  At  =  0.000125  is  chosen  for  .subsequent  calculations.  This  time  step  should  provide  adecpiale  temporal 
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resolution,  with  approximately  300  time  steps  per  cycle  at  the  highest  frequency.  No  siihiterations  ^r<- 
found  to  be  necessary  for  this  small  time  step. 


5.2  Low  Reynolds  Number  Case 

C'alculations  are  also  made  for  a  low  Reynolds  number,  Re  =  50,000.  This  provides  several  dist  inct 
advantages.  First,  the  calculations  are  in  the  Reynolds  number  range  of  the  experiments  of  Gad-el-Hak  [1] 
and  Lowson  [6].  Better  comparison  of  the  computed  results  with  experiment  is  then  po.s.sible.  Secondly, 
this  lower  value  of  Reynolds  number  also  eliminates  questions  involving  transition  and  turbulence  that 
ari.se  at  the  higher  Reynolds  number.  Finally,  the  lower  value  of  Reynolds  number  results  in  improved 
simulations  for  a  given  grid  size. 

The  delta  wing  geometry  itself  is  modelled  as  a  flat  plate.  This  is  done  to  simplify  the  grid  structure 
in  the  leading  edge  region  of  the  delta  wing.  It  is  felt  that  the  underside  geometry  has  no  significant  effect 
on  the  flow  physics  being  investigated. 

C'alculations  for  the  lower  value  of  Reynolds  number  are  carried  out  on  a  96  x  151  x  171  grid.  This 
grid  contains  54  points  in  the  axial  direction  and  80  points  in  the  spanwise  direction  on  the  delta  uiiq^ 
with  133  points  located  above  the  delta  wing.  The  minimum  spacing  at  the  wall  is  A:  =  O.UUOl  and 
the  spacing  along  the  leading  edge  of  the  delta  wing  varies  from  Aj/  =  0.0005  at  the  trailing  edge  to 
Ay  =:  3.165  X  10~®  at  the  apex.  The  spacing  in  the  axial  direction  at  the  apex  is  Ax  =  0.01  and  at 
the  trailing  edge.  A*  =  0.003.  For  the  predominant  portion  of  the  wing  there  is  a  constant  axial  sparing 
of  Ax  =  0.025.  This  grid  provides  significant  improvements  in  spatial  resolution  in  all  three  coordinate 
directions. 

The  unsteadiness  observed  at  high  values  of  Reynolds  number  remains  at  Re  =  50,000.  The  unsteady 
lielmvior  of  the  flow  field  for  this  case  is  represented  in  figures  7  and  8.  Each  figure  covers  a  time  period 
corresponding  to  approximately  one  cycle  of  the  unsteady  behavior  based  on  the  frequency  at  the  axial 
location,  x  =  0.7. 

In  figure  7  contours  of  the  a;-component  of  vorticity  are  shown  for  time  f  =  0.4375  to  i  =  0.4875  on  cross 
planes  located  from  0.3  to  0.85.  Positive  valu-  s  of  vorticity  are  represented  by  the  colors  yellow,  red.  and 
magenta,  with  magenta  representing  the  largest  values.  Negative  values  of  vorticity  are  repre.sented  by  (he 
colors  blue,  cyan,  and  green  with  blue  representing  the  smallest  values.  Clearly  visible  in  the  picture  is  lie 
primary  vortex  over  the  delta  wing.  Also  visible  are  a  series  of  di.stinct  small-scale  vortical  structure.s  hi 
the  shear  layer  which  emanates  from  the  leading  edge  of  the  delta  wing.  These  three-dimen.sional.  vortical 
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structures  are  first  evident  in  the  flow  at  an  axial  location  x  «  0.3.  The  small-scale  vortices  consist  of  a 
coherent  structure  that  forms  first  nearer  the  apex  of  the  delta  wing.  The  structure  does  not  extend  iiiitiall.v 
as  a  single  vortex  over  the  whole  length  of  the  leading  edge  of  the  delta  wing.  As  the  vortex  sheds  nearer 
the  apex  of  the  delta  wing  new  parts  of  the  vortical  structure  are  forming  further  downstream.  When 
tlie  vortex  sheds,  it  is  converted  around  the  primary  vortex  and  its  strength  is  dissipated.  This  pron'.s.s 
proceeds  much  more  rapidly  nearer  the  apex  of  the  delta  wing  where  the  scale  of  tlie  primary  vortex  is 
smaller.  As  the  previous  vortex  is  shed  a  new  vortex  is  formed  to  take  its  place.  The  temporal  evolution 
of  this  process  is  seen  in  figure  7.  These  small-scale  vortices  are  believed  to  l>e  the  source  of  the  striations 
ol).served  by  Gad-el-Hak,  figure  1,  and  also  observed  by  Lowson  [6]  in  his  smoke  flow  visualizations. 

Further  insight  into  the  ujisteady  behavior  may  be  obtained  by  looking  at  the  evolution  of  the  vortical 
structure  for  a  typical  cross  plane,  x  =  0.7,  figure  8.  In  this  figure  contours  of  the  axial  conipoiieiit  of 
vorticity  are  plotted  with  blue  contours  being  positive  and  red  contours  being  negative.  At  t  ime  i  =  t)..12.'). 
the  cycle  begins  with  a  small  scale  vortex  starting  to  form.  As  time  progre.sses  from  i  =  0.437.')  to  t  =  (l.dfl'i.*). 
this  vortex  grows  in  size  and  strength.  Subisequently,  for  t  ~  0.462.')  to  t  =  0.487.')  the  small-scale  vorti’X  is 
shell,  The  vortex  convects  around  the  primary  vortex  and  is  dissipated.  During  this  prores.s  of  slieiir-lnyer 
roll  u|)  and  vortex  shedding  an  interaction  can  be  observed  with  the  secondary  region  of  vorticity.  \'orticil> 
of  the  opposite  sign  is  shredded  from  this  region  and  is  dissipated  along  with  the  small-scale  vorlex  a.s  it  i.s 
convecteil  along  the  primary  vortex.  Even  though  contours  of  vorticity  may  not  be  compared  direct  ly  to 
smoke  flow  visualizations,  these  figures  appear  to  agree  qualitatively  with  those  of  Low.son.  figure  2. 

A  time  history  of  the  pressure  at  points  in  the  flow  field  where  the  passing  of  the  small-scalt‘  vorti'es 
may  be  monitored  is  made  at  several  axial  locations.  A  typical  time  history  for  the  axial  location  .r  =  (1.7 
is  given  in  figure  9.  A  dominant  frequency  is  clearly  seen  even  though  the  existence  of  other  harmoii'c.s 
is  evident.  This  dominant  frequency  is  found  to  correspond  to  the  shedding  of  the  small-scale  vortical 
structures.  While  the  time  histories  are  not  sufficiently  long  to  allow  for  an  accurate  frequency  domain 
analysis,  an  average  dominant  frequency  at  each  location  can  be  found,  and  the  results  are  shown  in 
figure  10.  The  average  frequencies  are  seen  to  scale  almost  linearly  with  axial  distance. 

The  effect  of  the  unsteady  behavior  on  the  .surface  pre.ssure  is  shown  for  axial  location  =  0.7,  figure  11. 
The  shedding  of  the  small-scale  vortices  contributes  to  only  a  small  temporal  variation  of  the  surface 
pressure  outboard  of  the  core  of  the  primary  vortex.  This  contrasts  with  the  high  Reynolds  number 
case  [8]  where  a  marked  temporal  effect  on  the  surface  pressure  outboard  of  the  core  of  the  primary  vortex 
was  seen.  For  the  high  Reynolds  number  the  small-.scale  vortices  form  and  are  shed  from  the  leading  <'d,4e 
of  the  delta  wing.  In  the  low  Reynolds  number  case  the  .small-.scale  vortices  form  further  away  from  the 


Ii'adiiig  txlgt*  of  the  delta  \vii)g,  thus  their  interaction  with  the  secondary  separation  region  and  with  the 
flow  near  the  surface  of  the  wing  is  reduced.  Time  Iii.stories  of  the  lift  and  drag  coefficients  also  show  only 
a  slight  temporal  variation,  figures  12  and  13.  The  lift  coefficient  varies  by  appro.ximately  2.5%  ainl  the 
drag  coefficient  varies  by  appimimately  1.8'%.  This  temporal  variation  occurs  at  a  much  lower  frei|uenc\ 
than  the  shedding  frequency  of  the  small-.scale  vortices. 

A  typical  simulated  surface  oil  flow  pattern  for  the  delta  wing  is  given  in  figure  14.  A  rather  e.\leii.'>ive 
region  of  reversed  flow  is  seen  outboard  of  the  secondary  separation  line.  This  reversed  flow  region  i*.\teiids 
to  nearly  the  midchord  of  the  delta  wing  in  an  area  located  between  the  secondary  and  tertiary  lines  of 
st'paration.  The  threi'-dimcnsional  extent  of  this  reversed  flow  region  is  .seen  in  figure  15  where  the  siirfr.-'- 
of  zero  .r-coinponent  of  velocity  is  plotted.  From  this  figure  it  is  seen  that  the  largest  region  of  revers.^d 
flow  is  located  between  the  75'%  chord  location  and  the  trailing  edge.  The  reversed  flow  region  e.Nli'iuliiig 
forward  to  approximately  midchord  is  only  a  small  tertiary  separation  region  near  the  surface  of  thi-  ih-lla 
wing. 
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6.  Discussion 


The  current  calculations  show  some  similarities  to  the  flow  visualizations  of  both  Ciad-el-tliik  and 
Blackwelder  [1]  and  Lowson  [(5].  Direct  comparison  of  the  computations  with  experiments  remains  tlillicnlt. 
hoWev(‘r.  and  only  qualitative  observations  are  possible.  The  majority  of  the  results  of  Lowson,  for  exainpb', 
are  subject  to  a  50  Hz  wind-tunnel  background  excitation  which  is  not  modelled  in  the  compulations. 
Furthermore,  comparisons  of  vorticity  contours  with  both  smoke  and  dye  flow  visualizations  are  iiol  iwact 
and  can  sometimes  be  misleading  [24],  Finally,  the  experiments  and  the  computations  are  all  done  for 
different  sweep  angles.  With  these  cautions  noted  the  following  observations  can  be  niaile. 

Hot  h  Gad-el-Hak  and  Lowson  view  the  instability  as  a  series  of  vortices  formed  and  shed  from  I  he  leading 
edge  of  the  delta  wing  at  a  .single  frequency.  Gad-el-Hak  gives  an  empirical  relation  for  this  frei|iieiic\ : 


while  Lowson  gives  the  relation: 


St  = 


St  = 


1625 

2577 

v/j^ 


(ti.t) 


(62) 


This  de,scription  of  the  instability  is  not  completely  consistent  with  what  is  ob.served  in  the  calculations, 
however.  Here  the  small-scale  vortices  do  not  form  at  the  leading  edge  but  further  along  the  vortex  sheet. 
Furthermore,  the  vortices  do  not  roll  up  as  a  single  unit  along  the  entire  length  of  the  lending  edge,  but 
rather  form  first  nearer  the  apex  of  the  delta  wing  and  then  subsequently  farther  back  along  the  h'nding 
edge  of  the  wing.  In  the  case  of  the  experiments  of  Lowson  this  difference  in  character  may  be  attribul  nble 
to  the  50  Hz  wind-tunnel  excitation  of  the  flow  field  forcing  the  vortices  to  form  at  the  leading  edgi*. 

Hatlier  than  a  single  frequency,  an  almost  linear  variatiotj  of  the  frequency  of  formation  of  the  simdl- 
scale  vortices  is  observed,  with  higher  frequencies  nearer  the  apex  of  the  delta  wing.  For  lit  =  5().l)(l(), 
the  empirical  formula  of  Gad-el-Hak, equation  (6.1)  predicts  a  frequency.  St  =  7.2  while  Lowson 's  relation, 
equation  (6.2)  predicts  a  frequency,  St  =  11.5.  These  values  are  of  the  same  order  as  the  fr(>qnencies 
computed  in  the  aft  portion  of  the  delta  wing.  This  gives  confidence  that  the  compuleil  fre<|uencies  are 
not  out  of  line  witii  the  actual  physical  situation. 

Further  insight  into  the  shear-layer  instability  is  obtained  by  considering  an  iiiviscid.  linear  slabibiv 
analysis  for  a  two-dimensional  shear  layer  between  co-flowing  streams.  Monkewitz  and  Htii'iie  [2.5]  have 
consiilrred  analytically  the  cn.se  where  the  mean  velocity  profile  is  described  by: 


\vli(»re  A  =  (Af’/Sr)  is  a  measure  of  the  velocity  difference  across  the  layer  with  U  the  averag(>  velocit\  of 
the  two  streams  and  1/  =  Ay* where  y*  is  the  dimensional  length  and  is  the  vorticity  thickness.  For 
A  =  1  tlie  iiiaximnm  amplified  frequency  from  the  linear  spatial  .stability  analysis  was  found  (o  he  u.  =  (I  ’/’I 
where: 


_  K-  (2a-/) 
4  f 


((i  I) 


'I'his  value  does  not  depend  strongly  on  A. 

In  an  iinlividual  cro.s.s-plane  at  a  given  time,  the  flow  which  emanates  from  I  lie  leading  edge  of  a  delta 
wing  may  be  viewed  as  a  skewed  shear  layer.  This  shear  layer  may  be  separated  into  two  com|ioni‘nts. 
an  axial  component  in  the  j-dircction  and  a  shear  or  cros.s-flow  component.  To  a  first  approximation  'he 
stability  of  the  shear  component  is  assumed  to  be  the  dominant  factor  in  the  stabilitv  of  the  shear  laver.  If 
the  shedding  freiiaeiu ies  at  the  cross-planes  previously  considered  are  nondimensionalized  as  in  et|nalion 
(().  I),  values  ofu;  are  obtained  in  the  range  u,'  «  0.13-0.32  with  the  specific  value  ofu;  dependent  on  the  local 
sliear-lujer  properties  at  the  location  cho.sen  for  nondimeusionalization.  These  values  ofu;  are  of  the  sann> 
order  as  the  maximum  amplified  frequency  u;  =  0.21  found  in  the  spatial  stability  analysis.  Furlhermoi-e. 
a  value  of  iv  as  0.21  always  occurs  at  the  approximate  location  where  the  small  scale  vortices  are  forming 
in  the  shear  layer.  This  simplified  analysis  indicates  that  the  computed  flow  field  nn.sleadine.ss  results  I'roni 
a  shear-layer  instability.  Furthermore,  the  observed  frequencies  are  con.si.stent  with  a  two-dimensioir.l. 
spatial,  linear  stability  analysis  for  the.se  types  of  .shear-layer  flows. 

Pairing  of  the  small-scale  vortical  structures  has  been  observed  in  both  the  experiments  of  fJnd-el-llak 
and  Lowson.  tlad-el-IIak  sttggests  that  the  primary  vortex  over  the  delta  wing  originates  n.s  a  .sei-i'.-  of 
pairings  of  the  smaller  scale  vc.tices  shed  from  the  leading  edge  of  the  delta  wing.  Pairing  of  the  sim.ll 
vortices  is  not  observed  in  the  computations.  The  reason  for  this  discrepancy  is  not  fully  understood.  The 
shear  layer  which  emanates  from  the  wing  leading  edge  is  expected  to  be  highly  susceptible  to  modification 
line  to  external  excitation,  intentional  or  unintentional.  This  excitation  of  the  shear  layer  could  pro."ole  the 
pairing  prore.ss  and  is  currently  not  modelled  in  the  computations.  In  order  to  obtain  a  rigorous  comparison 
of  computations  with  experiment,  a  known  forcing  frequency  might  be  required  in  both  experiment  and 
computation.  The  effect  on  the  pairing  process  due  to  numerica!  dissipation  in  the  calculations  is  also  not 
fully  understood,  and  re(|uires  further  investigation. 

It  should  also  be  noted  that  both  experiments  [26]  and  computations  [27]  for  flow  over  sleiulei  boiliesof 
revolution  at  large  incidence  exhibit  a  .similar  type  of  .shear-layer  instability.  Small  scale,  three-dinn’ii'-ion.'il 
vortices  are  oh, served  to  move  along  the  pumar\  .surface  of  cro.s.sflow  separation  emanating  from  the  bodv. 
'File  high  frequency  phenomenon  observed  is  distinct  from  the  \on  Karman-t\pe  siiediliiig  that  ni  rurs  •! 
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very  high  inridoiicf  for  s  'ncli'r  bodies  of  revolution.  The  experiments  [26]  and  computations  [27]  compare 
only  (|ualilntively.  but  the  numerical  calculations  capture  the  pliysical  instability. 

Finally,  computations  for  the  same  delta  wing  and  flow  properties  have  been  performed  In  others  using 
several  different  numerical  techniciues  [23.28.29].  These  computations  all  show  good  agreement  with  the 
experimental  measurements.  No  unsteady  flow  behavior  similar  to  the  instability  previonslj  <li.scn.vsed  i.s 
reported  in  the.se  references,  however.  Murman  [30]  has  observed  an  unsteaily  flow  behavior  that  results  in 
a  lack  of  convergence  of  the  re.sidual  to  steady  .state  for  Euler  calculations  over  della  wings.  The  .source  i.f 
the  lack  of  convergence  is  traced  to  the  existence  of  unsteady,  small-scale  vortices  in  the  shear  la\er  from 
the  wing  leading  edge.  .Since  the  computations  are  not  fully  time  accurate,  these  observations  are  onl\ 
(|iiaiitative.  A  full  understanding  of  the  effect  of  numerical  technique,  grid  topology,  and  griti  re.solulii>n 
on  the  compulation  of  this  tyjie  of  shear-layer  instability  will  require  further  investigations. 


7.  Conclusions 


A  iitnnerical  investigation  of  the  unsteady  vortex  structure  over  a  76®  sweep  delta  wing  at  20.5®  angle 
of  attack  has  been  carried  out.  Calculations  are  performed  for  a  Mach  number  0.2,  and  a  Heynolcl.s  imiiilx'r 
50. 000.  The  numerical  calculations  show  that  the  shear  layer  which  emanates  from  the  leading  e<|ge  of 
the  delta  wing  is  subject  to  a  Kelvin-Helmholtz  type  instability.  .Small-scale,  tlnee-dimen-sional  voriic.d 
.structures  are  observed  in  the  shear  layer.  The  small-scale  vortices  consist  ol  a  coherent  struct  m  e  tb;)' 
forms  first  nearer  the  apex  of  the  delta  wing  and  subsequently  farther  aft  along  the  leading  edge.  These 
vortices  are  shed  and  convected  around  the  primary  vortex  with  their  strength  being  di.ssipal(‘d.  The 
fre«iuency  of  shedding  of  these  small-scale  .structures  is  seen  to  vary  nearly  linearly  with  axial  location. 
'I'he  ob.sorved  frequencies  are  consistent  with  the  maximum  amplified  frequency  found  from  an  inviscid. 
linear,  spatial  stability  analy.sis  for  a  two-dimen.sional.  ro.ss-flow  shear  layer.  Ihe  small-.scale  vortical 
striictures  have  been  shown  to  be  related  qualitatively  to  those  observed  by  Gad-el-Ilak  [1]  ami  Low.son  [6] 
in  their  flow  visualization  experiments,  though  .some  differences  in  the  specific  character  of  the  computed 
and  experimental  structures  have  been  noted. 
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8.  Nomenclature 


f-v 

= 

pressure  coelficient 

t', 

= 

total  energy 

f 

= 

frecpiency 

■J 

transfonnaliou  Jacobian 

L 

= 

root,  chord  lengtii 

= 

‘‘reestreani  Mach  number 

l> 

= 

y.ressiire 

Pv 

= 

Prandtl  number 

fir 

= 

Reynolds  number 

.S7 

= 

Strouhal  number  St  =  fL/ti^ 

i 

= 

nondimensional  time  t  =  hi^/f. 

ti.  c.  «• 

= 

velocity  components  in  x,y  and  ; 

iV.  ■: 

physical  coordinates 

.</< 

= 

location  of  leading  edge 

CnX 

= 

computational  coordinates 

K- 

= 

vorticity  thickness 

p 

= 

density 

p 

= 

viscoisity  coefficient 

u> 

= 

nondimensional  frequency 

ri.j 

stre.ss  tensor 
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Figtirp  4:  Grid  Structure 
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t=0,475 


Figure  7:  Three-Dimensional,  Unsteady  Vortex  Structure  x  =  0.3  —  0.85, 
Contours  of  Axial  Component  of  Vorticity 
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